Exome Sequencing and Complex Disease: Practical Aspects of Rare Variant Association Studies

Genetic association and linkage studies can provide insights into complex disease biology, guiding the development of new diagnostic and therapeutic strategies. Over the past decade, genetic association studies have largely focused on common, easy to measure genetic variants shared between many individuals. These common variants typically have subtle functional consequence and translating the resulting association signals into biological insights can be challenging. In the last few years, exome sequencing has emerged as a cost-effective strategy for extending these studies to include rare coding variants, which often have more marked functional consequences. Here, we provide practical guidance in the design and analysis of complex trait association studies focused on rare, coding variants

Genome-wide association and linkage analyses of hemostatic factors and hematological phenotypes in the Framingham Heart Study

<p>Abstract</p> <p>Background</p> <p>Increased circulating levels of hemostatic factors as well as anemia have been associated with increased risk of cardiovascular disease (CVD). Known associations between hemostatic factors and sequence variants at genes encoding these factors explain only a small proportion of total phenotypic variation. We sought to confirm known putative loci and identify novel loci that may influence either trait in genome-wide association and linkage analyses using the Affymetrix GeneChip 100K single nucleotide polymorphism (SNP) set.</p> <p>Methods</p> <p>Plasma levels of circulating hemostatic factors (fibrinogen, factor VII, plasminogen activator inhibitor-1, von Willebrand factor, tissue plasminogen activator, D-dimer) and hematological phenotypes (platelet aggregation, viscosity, hemoglobin, red blood cell count, mean corpuscular volume, mean corpuscular hemoglobin concentration) were obtained in approximately 1000 Framingham Heart Study (FHS) participants from 310 families. Population-based association analyses using the generalized estimating equations (GEE), family-based association test (FBAT), and multipoint variance components linkage analyses were performed on the multivariable adjusted residuals of hemostatic and hematological phenotypes.</p> <p>Results</p> <p>In association analysis, the lowest GEE p-value for hemostatic factors was p = 4.5*10^-16for factor VII at SNP rs561241, a variant located near the <it>F7 </it>gene and in complete linkage disequilibrium (LD) (r²= 1) with the Arg353Gln <it>F7 </it>SNP previously shown to account for 9% of total phenotypic variance. The lowest GEE p-value for hematological phenotypes was 7*10^-8at SNP rs2412522 on chromosome 4 for mean corpuscular hemoglobin concentration. We presented top 25 most significant GEE results with p-values in the range of 10^-6to 10^-5for hemostatic or hematological phenotypes. In relating 100K SNPs to known candidate genes, we identified two SNPs (rs1582055, rs4897475) in erythrocyte membrane protein band 4.1-like 2 (<it>EPB41L2</it>) associated with hematological phenotypes (GEE p < 10^-3). In linkage analyses, the highest linkage LOD score for hemostatic factors was 3.3 for factor VII on chromosome 10 around 15 Mb, and for hematological phenotypes, LOD 3.4 for hemoglobin on chromosome 4 around 55 Mb. All GEE and FBAT association and variance components linkage results can be found at <url>http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007</url></p> <p>Conclusion</p> <p>Using genome-wide association methodology, we have successfully identified a SNP in complete LD with a sequence variant previously shown to be strongly associated with factor VII, providing proof of principle for this approach. Further study of additional strongly associated SNPs and linked regions may identify novel variants that influence the inter-individual variability in hemostatic factors and hematological phenotypes.</p

Genome-Wide Association with Select Biomarker Traits in the Framingham Heart Study

BACKGROUND: Systemic biomarkers provide insights into disease pathogenesis, diagnosis, and risk stratification. Many systemic biomarker concentrations are heritable phenotypes. Genome-wide association studies (GWAS) provide mechanisms to investigate the genetic contributions to biomarker variability unconstrained by current knowledge of physiological relations. METHODS: We examined the association of Affymetrix 100K GeneChip single nucleotide polymorphisms (SNPs) to 22 systemic biomarker concentrations in 4 biological domains: inflammation/oxidative stress; natriuretic peptides; liver function; and vitamins. Related members of the Framingham Offspring cohort (n = 1012; mean age 59 ± 10 years, 51% women) had both phenotype and genotype data (minimum-maximum per phenotype n = 507–1008). We used Generalized Estimating Equations (GEE), Family Based Association Tests (FBAT) and variance components linkage to relate SNPs to multivariable-adjusted biomarker residuals. Autosomal SNPs (n = 70,987) meeting the following criteria were studied: minor allele frequency ≥ 10%, call rate ≥ 80% and Hardy-Weinberg equilibrium p ≥ 0.001. RESULTS: With GEE, 58 SNPs had p < 10-6: the top SNPs were rs2494250 (p = 1.00*10-14) and rs4128725 (p = 3.68*10-12) for monocyte chemoattractant protein-1 (MCP1), and rs2794520 (p = 2.83*10-8) and rs2808629 (p = 3.19*10-8) for C-reactive protein (CRP) averaged from 3 examinations (over about 20 years). With FBAT, 11 SNPs had p < 10-6: the top SNPs were the same for MCP1 (rs4128725, p = 3.28*10-8, and rs2494250, p = 3.55*10-8), and also included B-type natriuretic peptide (rs437021, p = 1.01*10-6) and Vitamin K percent undercarboxylated osteocalcin (rs2052028, p = 1.07*10-6). The peak LOD (logarithm of the odds) scores were for MCP1 (4.38, chromosome 1) and CRP (3.28, chromosome 1; previously described) concentrations; of note the 1.5 support interval included the MCP1 and CRP SNPs reported above (GEE model). Previous candidate SNP associations with circulating CRP concentrations were replicated at p < 0.05; the SNPs rs2794520 and rs2808629 are in linkage disequilibrium with previously reported SNPs. GEE, FBAT and linkage results are posted at . CONCLUSION: The Framingham GWAS represents a resource to describe potentially novel genetic influences on systemic biomarker variability. The newly described associations will need to be replicated in other studies.National Heart, Lung, and Blood Institute's Framingham Heart Study (N01-HC25195); National Institutes of Health National Center for Research Resources Shared Instrumentation grant (1S10RR163736-01A1); National Institutes of Health (HL064753, HL076784, AG028321, HL71039, 2 K24HL04334, 1K23 HL083102); Doris Duke Charitable Foundation; American Diabetes Association Career Developement Award; National Center for Research Resources (GCRC M01-RR01066); US Department of Agriculture Agricultural Research Service (58-1950-001, 58-1950-401); National Institute of Aging (AG14759

Recurring exon deletions in the haptoglobin (HP) gene associate with lower blood cholesterol levels

Two exons of the human haptoglobin (HP) gene exhibit copy number variation that affects HP multimerization and underlies one of the first protein polymorphisms identified in humans. The evolutionary origins and medical significance of this polymorphism have been uncertain. Here we show that this variation has likely arisen from the recurring reversion of an ancient hominin-specific duplication of these exons. Though this polymorphism has been largely invisible to genome-wide genetic studies to date, we describe a way to analyze it by imputation from SNP haplotypes and find among 22,288 individuals that these HP exonic deletions associate with reduced LDL and total cholesterol levels. We show that these deletions, and a SNP that affects HP expression, are the likely drivers of the strong but complex association of cholesterol levels to SNPs near HP. Recurring exonic deletions in the haptoglobin gene likely enhance human health by lowering cholesterol levels in the blood

Clinical and Genetic Association of Serum Paraoxonase and Arylesterase Activities With Cardiovascular Risk

Objective—Diminished serum paraoxonase and arylesterase activities (measures of paraoxonase-1 [PON-1] function) in humans have been linked to heightened systemic oxidative stress and atherosclerosis risk. The clinical prognostic use of measuring distinct PON-1 activities has not been established, and the genetic determinants of PON-1 activities are not known. Methods and Results—We established analytically robust high-throughput assays for serum paraoxonase and arylesterase activities and measured these in 3668 stable subjects undergoing elective coronary angiography without acute coronary syndrome and were prospectively followed for major adverse cardiovascular events (MACE= death, myocardial infarction, stroke) over 3 years. Low serum arylesterase and paraoxonase activities were both associated with increased risk for MACE, with arylesterase activity showing greatest prognostic value (quartile 4 versus quartile 1; hazard ratio 2.63; 95% CI, 1.97–3.50; P\u3c0.01). Arylesterase remained significant after adjusting for traditional risk factors, C-reactive protein, and creatinine clearance (hazard ratio, 2.20; 95% CI, 1.60–3.02; P\u3c0.01), predicted future development of MACE in both primary and secondary prevention populations, and reclassified risk categories incrementally to traditional clinical variables. A genome-wide association study identified distinct single nucleotide polymorphisms within the PON-1 gene that were highly significantly associated with serum paraoxonase (1.18×10-303) or arylesterase (4.99×10−116) activity but these variants were not associated with either 3-year MACE risk in an angiographic cohort (n=2136) or history of either coronary artery disease or myocardial infarction in the Coronary Artery Disease Genome-Wide Replication and Meta-Analysis consortium (n≈80 000 subjects). Conclusion—Diminished serum arylesterase activity, but not the genetic determinants of PON-1 functional measures, provides incremental prognostic value and clinical reclassification of stable subjects at risk of developing MACE

Clinical and Genetic Association of Serum Paraoxonase and Arylesterase Activities With Cardiovascular Risk

Objective—Diminished serum paraoxonase and arylesterase activities (measures of paraoxonase-1 [PON-1] function) in humans have been linked to heightened systemic oxidative stress and atherosclerosis risk. The clinical prognostic use of measuring distinct PON-1 activities has not been established, and the genetic determinants of PON-1 activities are not known. Methods and Results—We established analytically robust high-throughput assays for serum paraoxonase and arylesterase activities and measured these in 3668 stable subjects undergoing elective coronary angiography without acute coronary syndrome and were prospectively followed for major adverse cardiovascular events (MACE= death, myocardial infarction, stroke) over 3 years. Low serum arylesterase and paraoxonase activities were both associated with increased risk for MACE, with arylesterase activity showing greatest prognostic value (quartile 4 versus quartile 1; hazard ratio 2.63; 95% CI, 1.97–3.50; P\u3c0.01). Arylesterase remained significant after adjusting for traditional risk factors, C-reactive protein, and creatinine clearance (hazard ratio, 2.20; 95% CI, 1.60–3.02; P\u3c0.01), predicted future development of MACE in both primary and secondary prevention populations, and reclassified risk categories incrementally to traditional clinical variables. A genome-wide association study identified distinct single nucleotide polymorphisms within the PON-1 gene that were highly significantly associated with serum paraoxonase (1.18×10-303) or arylesterase (4.99×10−116) activity but these variants were not associated with either 3-year MACE risk in an angiographic cohort (n=2136) or history of either coronary artery disease or myocardial infarction in the Coronary Artery Disease Genome-Wide Replication and Meta-Analysis consortium (n≈80 000 subjects). Conclusion—Diminished serum arylesterase activity, but not the genetic determinants of PON-1 functional measures, provides incremental prognostic value and clinical reclassification of stable subjects at risk of developing MACE

Relationship between moderate-to-vigorous physical activity, abdominal fat and immunometabolic markers in postmenopausal women

AbstractObjectsTo assess the burden of levels of physical activity, non-esterified fatty acids (NEFA), triacylglycerol and abdominal fat on the immunometabolic profile of postmenopausal women.Study designForty-nine postmenopausal women [mean age 59.43 (standard deviation 5.61) years] who did not undertake regular physical exercise participated in this study. Body composition was assessed using dual-energy X-ray absorptiometry, and levels of NEFA, tumour necrosis factor-α, adiponectin, insulin and triacylglycerol were assessed using fasting blood samples. The level of physical activity was assessed using an accelerometer (Actigraph GTX3x), and reported as counts/min, time spent undertaking sedentary activities and time spent undertaking moderate-to-vigorous physical activity (MVPA). The following conditions were considered to be risk factors: (i) sedentary lifestyle (<150min of MVPA per week); (ii) high level (above median) of abdominal fat; and (iii) hypertriacylglycerolaemia (<150mg/dl of triacylglycerol).ResultsIn comparison with active women, sedentary women had higher levels of body fat (%) (p=0.041) and NEFA (p=0.064). Women with higher levels of abdominal fat had impaired insulin resistance (HOMA-IR) (p=0.016) and spent more time undertaking sedentary activities (p=0.043). Moreover, the women with two risk factors or more had high levels of NEFA and HOMA-IR (p<0.05), as well as an eight-fold higher risk of a high level of NEFA, independent of age (p<0.05). No significant relationship was found between levels of physical activity, abdominal fat, tumour necrosis factor-α and adiponectin (p>0.05).ConclusionPostmenopausal women with a combination of hypertriacylglycerolaemia, a high level of abdominal fat and a sedentary lifestyle are more likely to have metabolic disturbances

Fabrication of Tizanidine Loaded Patches Using Flaxseed Oil and Coriander Oil as a Penetration Enhancer for Transdermal Delivery

Transdermal drug delivery is important to maintain plasma drug concentrations for therapeutic efficacy. The current study reports the design, formulation, and evaluation of tizanidine transdermal patches formulated using chitosan and thiolated chitosan, ethyl cellulose (EC), polyvinylpyrrolidone (PVP), and Eudragit RL100 in different ratios. The tizanidine patches were formulated using flaxseed oil and coriander oil in the concentrations of 1% v/w, 2% v/w, 3% v/w, 4% v/w, 5% v/w, and 10% v/w. The patches were subjected to characterization of physicochemical property (thickness, weight uniformity, drug content, efficiency, percentage moisture uptake/loss), in vitro drug release and drug permeation, skin irritation, in vivo application, pharmacokinetics analysis, and stability studies. The results indicate that the interaction of thiolated chitosan with the negative charges of the skin opens the tight junctions of the skin, whereas flaxseed and coriander oils change the conformational domain of the skin. The novelty of this study is in the use of flaxseed and coriander oils as skin permeation enhancers for the formulation of tizanidine transdermal patches. The formulations follow non-Fickian drug release kinetics. The FTZNE23, FTZNE36 and FTZNE54, with 5% v/w flaxseed oil loaded formulations, exhibited higher flux through rabbit skin compared with FTZNE30, FTZNE35, FTZNE42, and FTZNE47, formulations loaded with 10% v/w coriander oil. The study concludes that flaxseed oil is a better choice for formulating tizanidine patches, offering optimal plasma concentration and therapeutic efficacy, and recommends the use of flaxseed and coriander oil based patches as a novel transdermal delivery system for tizanidine and related classes of drugs

Genome-wide association of echocardiographic dimensions, brachial artery endothelial function and treadmill exercise responses in the Framingham Heart Study

<p>Abstract</p> <p>Background</p> <p>Echocardiographic left ventricular (LV) measurements, exercise responses to standardized treadmill test (ETT) and brachial artery (BA) vascular function are heritable traits that are associated with cardiovascular disease risk. We conducted a genome-wide association study (GWAS) in the community-based Framingham Heart Study.</p> <p>Methods</p> <p>We estimated multivariable-adjusted residuals for quantitative echocardiography, ETT and BA function traits. Echocardiography residuals were averaged across 4 examinations and included LV mass, diastolic and systolic dimensions, wall thickness, fractional shortening, left atrial and aortic root size. ETT measures (single exam) included systolic blood pressure and heart rate responses during exercise stage 2, and at 3 minutes post-exercise. BA measures (single exam) included vessel diameter, flow-mediated dilation (FMD), and baseline and hyperemic flow responses. Generalized estimating equations (GEE), family-based association tests (FBAT) and variance-components linkage were used to relate multivariable-adjusted trait residuals to 70,987 SNPs (Human 100K GeneChip, Affymetrix) restricted to autosomal SNPs with minor allele frequency ≥0.10, genotype call rate ≥0.80, and Hardy-Weinberg equilibrium p ≥ 0.001.</p> <p>Results</p> <p>We summarize results from 17 traits in up to 1238 related middle-aged to elderly men and women. Results of all association and linkage analyses are web-posted at <url>http://ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?id=phs000007</url>. We confirmed modest-to-strong heritabilities (estimates 0.30–0.52) for several Echo, ETT and BA function traits. Overall, p < 10^-5in either GEE or FBAT models were observed for 21 SNPs (nine for echocardiography, eleven for ETT and one for BA function). The top SNPs associated were (GEE results): LV diastolic dimension, rs1379659 (<it>SLIT2</it>, p = 1.17*10^-7); LV systolic dimension, rs10504543 (<it>KCNB2</it>, p = 5.18*10^-6); LV mass, rs10498091 (p = 5.68*10^-6); Left atrial size, rs1935881 (<it>FAM5C</it>, p = 6.56*10^-6); exercise heart rate, rs6847149 (<it>NOLA1</it>, p = 2.74*10^-6); exercise systolic blood pressure, rs2553268 (<it>WRN</it>, p = 6.3*10^-6); BA baseline flow, rs3814219 (<it>OBFC1</it>, 9.48*10^-7), and FMD, rs4148686 (<it>CFTR</it>, p = 1.13*10^-5). Several SNPs are reasonable biological candidates, with some being related to multiple traits suggesting pleiotropy. The peak LOD score was for LV mass (4.38; chromosome 5); the 1.5 LOD support interval included <it>NRG2</it>.</p> <p>Conclusion</p> <p>In hypothesis-generating GWAS of echocardiography, ETT and BA vascular function in a moderate-sized community-based sample, we identified several SNPs that are candidates for replication attempts and we provide a web-based GWAS resource for the research community.</p